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ABSTRACT 

The analysis of how a stochastic background of gravitational radiation interacts with a 
spherical detector is given in detail, which leads to explicit expressions for the system 
response functions, as well as for the cross- correlation matrix of different readout channels. 
It is shown that distinctive features of GW induced random detector excitations, relative 
to locally generated noise, are in practice insufficient to separate the signal from the noise 
by means of a single sphere, if prior knowledge on the GW spectral density is nil. The 
situation significantly improves when such previous knowledge is available, due to the 
omnidirectionality and multimode capacities of a spherical G W antenna. 
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1 Introduction 



Stochastic backgrounds of Gravitational Waves (GWs) are amongst the most interesting 
sources to be detected by the upcoming new generation of antennas: such kind of grav- 
itational radiation will convey to us e.g. information on the structure of the Universe at 
its earliest evolutionary stages. This is quite unique to gravitational waves, due to their 
extremely weak coupling to other forms of matter and radiation, and will therefore greatly 
enhance our understanding of Cosmology as well as Fundamental Physics. Literature on 
the possible forms, origins and spectral shapes of various gravitational wave backgrounds 
is rich — see for example for a detailed review and further references on the subject. 

The stochastic nature of gravitational background signals, together with the same 
stochastic origin of the local noise generated in the detector itself, makes very difficult 
to tell local noise from actual signal. Therefore, detection strategies to extract that kind of 
gravitational signals from detector data are based on finding differential properties capable 
to distinguish between both processes. 

The standard procedure to identify a stochastic signal is to use two (or more) anten- 
nas, and cross-correlate their outputs. The concept is that, while local noises in separate 
detectors are uncorrelated, the stochastic signal is common to both. Signal to noise ratio 
therefore builds up on the basis of long-term average compensation of local noise fluctu- 
ations, and can be seen to develop as T 1//4 , where T is the integration time [0. Specific 
strategies for GW background detection have been considered for two bar detectors [[|, 
two interferometers Q , a bar and an interferometer || and two spherical detectors || . 

A spherical GW detector is a multimode device which, in a number of senses, behaves 
like an array of bars ]7[ 0]. This suggests that one might try and take advantage 
of such multimode capacity to perform a search for a GW background by suitably cross- 
correlating the different detector readouts. A procedure like this would therefore enable the 
determination of the GW background spectral density with a single detector — a significant 
advantage to rid the usual method (sketched above) of inherent uncertainties bound to the 
fact that no two real detectors are exactly identical. 

One might a priori expect in this direction that specific signal correlation patterns 
happen between pairs of sensor readouts which be not shared by the local sources of noise, 
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thereby making possible to filter the latter out. This paper is concerned with the analysis 
of the possibilities offered by the sphere, as a multimode device, for the detection of a 
stochastic background of GWs, defined by a spectral density function. 

As we shall see, though, local noise and GW background random amplitudes actually 
show no distinctive correlation pattern in the output channels, except for the obvious 



fact that GWs are only seen in quadrupole (perhaps also monopole) antenna modes [10 



So, in the absence of some kind of previous information on the functional form of the 
GW spectrum, the above expectations fail. If, on the other hand, spectral information is 
available ahead of time (templates) then the sphere has a very good performance in its 
frequency band of sensitivity. This can improve efficiency up to an order of magnitude in 
energy compared with single readout GW detectors (bars and interferometers), due to the 
sphere being both omnidirectional and multimode, as we shall see. 

The outline of the article will be the following. In section ^| we set up definitions and 
notation for the stochastic background of GWs, and in section |3| we derive the sphere's 
response to this signal. In section |] we present the correlation functions between the 
sphere's responses induced by the gravitational signal at two points on its surface, where 
motion sensors will be attached. Section |5| is dedicated to calculate the same correlation 
function presented in section | but in this case induced by local noise, and discusses, in the 
light of these results, various possible strategies to retrieve information about the incoming 
gravitational signal spectrum with a single spherical detector. We close with a summary 
of conclusions in section 0, and also include a brief mathematical appendix. 



2 The stochastic signal 

We give in this section a few definitions and notation for a stochastic background of grav- 
itational radiation, and its characterization through the spectral density Sh(u). 

A background of gravitational waves can be expressed as a linear superposition of plane 
waves coming from all directions and with all possible polarizations, 

Mt,x) = T ^ f d 2 n e -^- n ^ r # £Ma,;n,v) Gj}(n,v) (i) 

J-oo Zn J Jo a 
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where f2 is a unit vector pointing to a generic location, if; is a rotation angle around 
this vector, 1ia(oj; ft, ijS) are the wave amplitudes and G^(ft,i^) the polarization matrices. 
General Relativity predicts only two transverse polarization modes. Nevertheless, another 
four degrees of freedom must generally be allowed for. We thus have altogether six in- 
dependent 3x3 symmetric polarization matrices. If we choose a right-handed triad of 
orthonormal vectors m, n and ft, where ft is aligned with the propagation direction, a 
suitable parametrization for the G-matrices is the following: 

g% (ft, i>) = \j\ tij , g§» (n, V) = \/| (3^ - <y 

Gff (SI, -0) = m 4 3 + , Gff (ft, tp) = nfrj + 

Gj* (ft, ip) = miUj + niirij , Gjj~(ft, ip) = miirij — riiUj (2) 

where we give to the index A in Gfj(ft, ip) the following meanings: A —M is the monopole 
component, A = Lq, L m , L n are longitudinal quadrupole components, and A = T + ,T X are the 
two transverse quadrupole components. These matrices satisfy normalization conditions 

Y.Gp% = 25 AB (3) 

ij 

If (9, 0, ip) are Euler angles relating the laboratory frame to the just described wave 
frame then 



m = (cos ip cos 9 cos — sin ip sin , cos ip cos 9 sin + sin ip cos , — cos ip sin 9) 
n = (— sin ip cos 6* cos — cos sin , — sin ip cos sin + cos ip cos , sin ip sin 0) 
ft = (sin cos 0, sin sin 0, cos 9) (A) 

where we can see that a choice for the angles 9, gives a direction in space and a choice for 
ip gives the definition of the two transverse polarization states for a particular direction. 

The " electric" components R 0i0 j(t) of the Riemann tensor at the center of mass of the 
sphere are 

RoiOj(t) = ^,oo(0, t) = - du u 2 e~^ J d 2 Q J # £ M"l ^ G£(fi, V), (5) 
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and these are the only components carrying information on an incoming GW JOJ ]. 

In this paper we shall be considering an unpolarized stochastic bath of GWs which is 
also isotropic and stationary. Its statistical properties are therefore encoded in its power 
spectrum function, defined by the following equation: 



(h* A (u,n,i/>)h A iu/,rt,1/)) = 6(u-<S) |^5 2 (n,n')^(^-^)} {<W'2^)} (6) 

where (— ) stands for ensemble average. There are alternative ways to characterize the 
frequency spectrum, for example, through an energy density per unit logarithmic inter- 
val of frequency, or in terms of a characteristic amplitude of the stochastic background. 
The definitions for such quantities, as well as their mutual relationships can be found in 
reference |]J. 



3 The sphere's response 

In this section we calculate the sphere's response to the background of gravitational waves 
given by equation (p]). We shall use henceforth the general formalism and notation of 



reference [10] 



We assume that the sphere's response u(x, t) to an incoming signal is the solution to 
the partial differential equation 

- { M V 2 u+(A + /i)V(V-u)} = f (x, t) (7) 

Cu 

with suitable initial and boundary conditions. Here, A and \i are the material's elastic 
Lame coefficients [12], and g is the sphere's density. In the rhs of equation ([?]), fj(x,i) = 



gc Roioj(t)xj is the density of gravitational wave tidal forces, ensuing from the geodesic 



deviation equation — see e.g. [13]. The latter term can be split up into its monopole and 



quadrupole components according to [[10] 

f(x,t) = f(°°)(x)/ )0 )(t) + £ f( 2m )(x)(/ 2m )(t) 
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with 



ff 0) (x) = qE^ x s , gM{t) = f E^ R^it) c 2 

ff m) (x) = e E% m) Xj , s (2m) (t) = ^ Ef m) R 0t0j (t) c 2 (m = -2, . . . , 2) (9) 

The definition for the matrices -E^ 00 ^ and E^ m ^ can be found in appendix A. The generic 
response function u(x, t) can be expressed as an orthogonal series expansion, which only 
involves quadrupole and monopole terms: 



u( x >0 = — u noo(x)^i° 0) (t) + £ — 

n=l ^0 n=l ^"2 



H U n2m(x)^2 m) (^) 



(10) 



(2m), 

m=-2 

where u n i m are spheroidal wavefunctions corresponding to the {21 + l)-degenerate eigen- 
frequency u nh i.e., 

£u nZm (x) = -u) 2 nl Q\i nlrn (x) (11) 

with £ the differential operator defined in equation (0); a n0 and a„2 are projection coeffi- 
cients — see [H for details on definitions and values of these quantities. 

Finally, <7^^ (£) are convolution products between the driving terms g^ lm \t) and the 
corresponding eigenmode oscillation function — in this case represented by an ideally non- 
dissipating, purely sinusoidal vibration: 

g$T\t)= fdt'g« m \t')sm{oo nl {t-t')} (1 = 0,2 ; m = -I, (12) 
Jo 

We shall later relax this ideal behavior hypothesis in order to cope with the stationary 
state the system will eventually reach under the action of long term random excitations. 

We next calculate the sphere's response to the background of gravitational waves given 
by equation (J^). The only quantities we have to evaluate are the functions g^^^t) using 
equations (0) and (|9|), with the components of the Riemann tensor given by equation fl5|). 
The following is readily found: 



9$T\t) = ~ r divu 2 T(t;u,u nl ) f rf 2 fi# ^h A (u;fl,i;) Gf^) E*f m) (13) 
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where 



T{t; u, uo nl ) = [ d£ e~ iuJt ' sin {u nl (t - t')} (14) 
Jo 

is a function strongly peaked at u — rku n i for large values of t. The sphere's response 
u(x, t) is then given by equation ([K]) with g^P (t) given by equation (pT3|) . 

In actual practice, the sphere's motions are sensed at a number of different positions 
on its surface, where motion sensors are attached. We consider devices which are only 
sensitive to radial displacements, and attached to the sphere's surface at positions 

x a = Rn a , a = 1, . . . , J (15) 

where R is the radius of the sphere, and n a a unit vector pointing outward at the a-th 
position. There will consequently be J readout channels with output displacements 

u a (t) = n a ■ u(x a , t) , a = l,...,J (16) 
If equation ([LC]) is now used, these are seen to be 



u 



GW 



Q-nO 



Y, — A n0 (R)Y 00 (n a )gl? 0) (t) 



+ E — A n2(R) 
n=l U n2 



+ 



E >Wn a )si 2 2 m) (t) 



1. 



J 



(17) 



where Y\ m are spherical harmonics, and A n i(R) is a radial function coefficient in the 
spheroidal wavefunction u nim (x); this one depends on whether we are dealing with a solid 
or a hollow sphere — see JT0| 



or 



|15|1 for each case, respectively. 



Because the time dependent terms gn^it) in (|T7| ) are random in nature, the readouts 
u a (t) are consequently random, too. In addition, the linear character of the relationship 
between both ensures that certain properties of the driving terms, such as e.g. gaussianity, 
directly carry over to the outputs. We now propose to investigate the relevant statistical 
properties of the latter. 
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4 Cross-correlation functions 



As already stated in the Introduction, an essential ingredient in our analysis are the cross- 
correlations between different output channels. In order to make them useful, however, 
we first need to establish their theoretical structure, i.e., to assess their behavior under 
ideal conditions of infinite integration times, and relate them to the unique characteristic 
function Sh(u) — the signal spectral density. 

Correlations between the outputs at the readout channels are naturally defined by 

R ab {t,r) = {u* a {t)u b {t + r)) (18) 

where (— ) stands for ensemble average again. We must now introduce the function u^ w (t), 
given by equation (|T7]), into equation (|T5J) to calculate the correlation functions we are 
looking for. When averages are taken, equation (|6|) for the ensemble average of the wave 
amplitudes comes into play. Dirac ^-functions appearing in that equation allow us to easily 
perform the integrals in f2', if)' and u 1 . The remaining integrals in f2, ip can then be done 
after evaluating 

72 r 



G m = J g / £ e^VOeftn,^) (19) 
which yields the following result: 



Gijki — - Sij5ki + 2N 



1 1 
15 10 



(20) 



monopole quadrupole 

The index A in equation ( |T9"D runs over the polarization modes of the incoming gravita- 
tional signal. And the result given in equation fl2~0|) takes into account the monopole mode 
and iV quadrupole modes (note that N < 5, for example iV = 2 in General Relativity, and 
that all of them have the same contribution to G^ki). 

It is now readily seen that there is no cross-correlation between monopole and quadru- 
pole modes, and that the correlation functions between the outputs induced by the gravi- 
tational signal split up as 

nS v (t,T) = Rfr\t,T ) + RT w (t,T) (2i) 
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where R a ° h ^ GW and i?^ Giy are monopole and quadrupole terms, respectively, and have the 
form 



RT w (t,r) = 2 P (n tt -n 6 ) ± -noa n/0 A n0 (R)AMR) &M 

n,n'=l U n0 UJn>0 

d(2)gw ( . s 4iV p (n n v ^ a n2 a n '2 A n2 (R)A n , 2 (R) ( 2 ) , \ (r>r> , 

10 n,n'=l W n 2U} n '2 

with Pq the zero-th Legendre polynomial — which is identically equal to 1 for any value of 
its argument, but we keep it explicit for the sake of structure clarity — , P 2 is the second 
Legendre polynomial, and 



,,n poo Atjj 

fW(t, t)= —LU 4 Re [T*{t, to, co nl ) T(t + r, u, S h {u) , I = 0, 2 (23) 
Jo Zn 

Here we see how the spectral density Sh(oo) characterizing the incoming signal appears, 
and its relation to the correlation functions. 

In a stationary random process, the second order correlation functions are time-shift 
invariant (by definition), i.e., if we use the definition (18) then the function R should 



only depend on its second argument, r. Time-shift invariance has been assumed for the 
exciting GW bath, as we have characterized it by the spectral density Sh(w)- Nevertheless, 
the sphere's response is calculated with the series expansion (|T0|) in which a choice for the 
state of the sphere at t — (initial conditions) is assumed which, therefore, breaks that 
invariance. 



The correlation functions we are looking for can actually be obtained from (|2T| ) and ( p2|) 
by taking their limit when t approaches infinity, because in this limit the property of 
time shift invariance is recovered, as the system looses memory of any particular initial 
conditions in the remote past. 

Memory is however not lost in the kind of ideal, non-dissipative elastic body we have 
described in the preceding sections. Meaningful results can only be obtained if some kind 
of dissipation, no matter how large or small, is present in the system. This is therefore 
the appropriate place to introduce dissipation, and we shall do it in the standard way of 
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assuming an exponentially decaying amplitude for the oscillation eigenmodes, i.e., we make 
the replacement 



sin(c^r) — ► e 7 " iT sin(cj n/ r) (24) 

where 7„/ is proportional to the inverse of the decay time of the associated mode, which is 
in all cases of interest much larger than the period of oscillation of that mode, or 

Inl < u n i (25) 

If we express the limits of R^ GW (t, r) and R^ GW (t, r) by dropping the first argument 
then the following results are readily obtained: 

«S ,G "'W= 2 fl>(SVn»)£^ fl«V) e*" (26) 

RT W (t) = H p 2(fia . n h ) r p (27) 

lb J-oo Z7T Z 

where 

fl^w) = 2 KiA nl {R)f u 4 \L nl {uj)\ 2 (28) 

n 

and L n i(u) is a Lorentzian curve: 

LmH = — ^i— (29) 

Equations (^6|) and (B7|) have a foreseeable functional structure: the correlation func- 
tions are Fourier transforms of the spectral density of the driving stochastic signal times 
the sphere's transfer function RW(u)), and this is a characteristic sum of Lorentzian curves 
centered on the resonant frequencies, and with appropriate weights. Note also that different 
harmonics (n, n') appear to be uncorrelated. 
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4.1 Mode channels 

The algebraic structure of the correlation function matrix requires some further consider- 
ation. For P2(Q a - ^fe) is a symmetric matrix whose rank[| is at most 5, which means one 
can form 5 linear combinations 

y (m) (*) = E«£ m) «.(*) . m = -2,...,2 (30) 

a=l 

where the coefficient is the a-th component of the (normalized) eigenvector associated 
to the m-th non-null eigenvalue of P 2 (fl a - fib), Cm' sa y' which is always positive. 
If we now define a new correlation function matrix 



= (y (m> (t)y^\t + r)) , m,m' = -2,...,2 (31) 

we readily see that 

«(t) = ^C^w r ^ B?\u) l -S h {u)e^ (32) 
15 J-oo 2n 2 

Thus the 5 quadrupole channels y^ m \t) are uncorrelated with one another. They are 
therefore particularly well suited for improved estimation of the spectral density Sh{oo), 
since the latter is common to all of them. One may recall here that y^ m '{t) are actually 



mode channels for the well known TIGA JT6| and PHC |HJ sensor distributions, in which 

cases 



v { r ] oc Y 2m {tl a ) for TIGA and PHC (33) 

Note that the fact that the y {m) {t) are uncorrelated is independent of whether the 

motion sensors are resonant or not, so it applies to a system like the recently proposed dual 

sphere ||18| , too. Note also that sensor geometries such as TIGA and PHC, for which the 

y( m \t) are mode channels — i.e., quantitites whose frequency spectrum y^ m \uj) is directly 

proportional to that of the GW quadrupole amplitudes g^ 2m \uj) — , allow for a very clear 

1 Actually, rank{P;(f2 Q - fit)} = 21 + 1, provided no two of the J vectors fl a are parallel, and provided 
of course that J > 21 + 1 — see H for full details. 
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and direct physical interpretation of the correlation coefficients (|32|). They are therefore 



also preferred from the specific point of view of the signals considered in this paper. 

Summing up, the main result so far is that cross-correlations between different output 
channels of a spherical GW detector, whether direct sensor readouts or quadrupole (mode) 
channels, possess a specific pattern, which is given by equations fl27j) or (|32|), respectively 
The problem is of course that a random signal must be told from also random local noise; 
we thus review next the relevant properties of the latter. 

5 Local noise and detection strategies 

We study in this section the multipole characteristics of the correlation functions between 
the outputs at two points of the sphere induced by the local noise. 

The actual readouts of the spherical antenna are in fact superpositions of those induced 
by the GW and by the local noise, i.e., 

u a (t)=u G a w (t)+u L a N (t) (34) 
We shall generically consider that 

< N (t)= E qnim(t)A nl (R)Y lm (n a ) (35) 

n, l,m 

where q n im{t) are narrow band stochastic time series, with correlation times of order l/7 n j. 
Because these q n im(t) are the amplitudes of different normal modes of the solid, they are 
statistically uncorrelated, or 

(ClmW Qn'l'm'(t + r)) = b nn i b w 8 mm i ^(r) (36) 

with, 

/DO fll,] 1 
*=- |L n ,H| 2 -S LN {u)e^ (37) 
-oo Z7T I 

where L n i{ui) is given by (p9|) , and Sln(uj) is the input spectral density of local detector 
noise. We need not go into the details of the specific form of this function here, which is a 
superposition of thermal and transducer back-action noise — see flTBfl. 



12 



The system readout correlation functions of local noise, defined by 



R L ab N (r)^(u L a N *(t)u^(t + r)) (3* 



are therefore given by 



R L ab N (r) = E RT N (t) (39) 



where 



RT N (t) = ^ m a - tl b ) jT ^ £ \A nl (R)\ 2 \L nl (^ \ S LN (") (40) 



These expressions show that local noise correlations split up in a way similar to (21), 
except of course that we now have contributions from all multipoles rather than just 
monopole and quadrupole modes, characteristic of GW signals. In particular, we can 
easily assess how do quadrupole channels — as defined in ( |30"D — correlate with one another 
in a local noise dominated detector: 



2 

+ E ^ POL % E I4u(*)| 2 IM-)| 2 ~ S LN {<*) (41) 



n2io;)| 2 -^jvM e WT 



where 



r$L = E vjrH^ fl(*V n 6 ) , m = -2, . . . , 2 (42) 



EJ.6 



is generally a non-diagonal matrixF[ 



The important conclusion to be drawn from equations ( flO|) and (41) is that monopole 



and quadrupole correlations follow precisely the same algebraic pattern, whether they are 
2 Exceptionally, it is diagonal for certain values of I in certain sensor configurations — see a discussion 
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dominated by local noise or by a bath of random GWs. Extra terms, such as those in the 
last line of (|4lD, do not really provide a distinctive feature of local noise, for they are only 
significant at frequencies different from those of the quadrupole modes, as determined by 
the presence of the frequency dependent coefficients L n i(u) in each case. 

This is a very unfortunate circumstance, indeed, for it in principle renders useless the 
cross-correlation algorithm between different detector channels as an efficient means of 
filtering out local noise. Are there expedient alternatives? 

5.1 Filtering strategies 

There are two interesting possibilities one can try to filter out local noise effects with a 
single multimode detector: 

i) Since different normal modes are statistically uncorrelated, one is tempted to consider 
e.g. the first and second quadrupole modes of a sphere || [l5|, |l^| as a pair of inde- 



pendent antennas, then use cross-correlation between them to filter local noise out. 

This however cannot possibly work because, even if local noise induced fluctuations 
are independent in either mode, so are also GW induced fluctuations. And, as we 
have just seen, there is no distinctive feature as to how noise gets added to the signal 
at different quadrupole modes. Actually, as seen in equations ( [31~1) and (f4lD , there is 



no way to tell which part in y^ m \t) is due to GWs and which is due to local noise. 

ii) The fundamental difference between a bath of random GWs and local noise is that, 
while the former can only affect quadrupole modes (or, at most, also monopole 
modes), local noise excites all modes, instead — not just quadrupole/monopole. 
Given that, in a realistic GW antenna, mode linewidths are extremely narrow, i.e., 
energy transfer between different modes is extremely slow, an attractive possibility 
to use a single sphere as a detector of a continuous flow of stochastic GWs would be 
to identify the effect of the latter as a temperature excess in the quadrupole modes, 
relative to other modes — dipole, octupole, etc. 

While this may look like a reasonable approach, it is in fact an impractical one, as we 
can appreciate by the following argument. The temperature of a mode is of course 
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a measure of the variance of its random oscillations; since these are a narrow band 
stochastic process, their correlation time, r corr , is quite long, hence a reliable estimate 
of the variance requires taking data during a suitable number of correlation times for 
averaging. If t b s is such time then the relative error in the temperature estimation 
is of the order of 



~ / 7 ~ corr (43) 
T y t bs 

Let us consider a few likely numbers: an SQL detector — such as described e.g. 
in |18| — should be capable of sensing energy innovations of magnitude kT / Q, where k 



is Boltzmann's constant, T is the temperature of the mode (e.g. the first quadrupole), 
and Q its mechanical quality factor. An optimized device should be able to reach 
T/Q ~ 10 _8 K _1 , for example with Q = 10 7 and T = 100 milli-Kelvin. A spherical 
antenna whose considered resonance happens at 1 kHz, say, thus has a correlation 
time of InQjioyi = 10 4 seconds. On the other hand, the precision required in the 
measurement of the mode temperature to match the SQL is ST/T = 10~ 7 for the 
assumed temperature of 0.1K. If we now make use of ( f43"D then the conclusion is 
that an integration time of about 3xl0 10 years is required to detect such a small 
temperature excess as the optimized antenna permits. . . 

It can be argued that stochastic background signals will actually be hotter than the 
just described precision limit. However nucleosynthesis constraints, for instance, set 
bounds on the spectral density of GWs to values below 3x 10 -24 Hz -1 / 2 at frequencies 
in the 1 kHz range [|IJ] . This translates into a GW bath temperature of fractions of a 
micro-Kelvin, for whose accurate estimation integration times of the order of several 
10 000 years would be required. The use of fast variance estimation techniques PU 



may improve on this by between one and two orders of magnitude, but even so one 
would still be in the range of 100 to 1000 years of integration time — an absurd figure. 

There is also little hope that GW backgrounds originating in astrophysical sources, 
such as e.g. supernovae [^TJ, have an effective temperature above 10 _7 K or so |1|]. 
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Altogether then the multimode capabilities of a spherical GW antenna do not seem to 
offer the practical possibility to make reliable, single detector measurements of the spectral 
density of a background of stochastic GWs. An exception to this however happens when 
some a priori information on the actual form of the spectral function is known ahead of 
time (template) ||22|| . If this is the case then optimum filter techniques can be used, of 



which a wealth is available in the literature 23 



Of course this is not specific to spherical detectors — in fact any single, unimodal 
GW antenna (such as bars or interferometers) can make advantageous use of the above 
techniques. Spheres however are more efficient in two senses: first they are omnidirectional, 
and second they can provide up to five independent estimates of the spectral density of 
the GW background, as follows from equation (|32|) : each quadrupole channel y^ m \t) can 
be used to produce one such independent estimate. On average, one gets a factor of 15/8 
because of omnidirectionality, times a factor of 5 for quadrupole channels, which make up 
for 75/8, one order of magnitude (in energy) better than bars and interferometers. Let us 
however stress that this obviously applies only within the frequency band of sensitivity of 



each detector; in the case of a dual sphere, even this is quite competitive, see [18 



6 Conclusion 

In this paper we have investigated in detail how a spherical detector interacts with a 
background of gravitational waves, characterized by a frequency dependent spectral energy 
density function, Sh(oo). In particular, we have found the generic response functions for 
a device with an arbitrary number of motion sensors, as well as their correlation matrix, 
both of the system readouts and for the five quadrupole channels — mode channels for 
certain specific sensor layout geometries. 

We have discovered that there is no distinctive pattern in that correlation matrix, 
relative to the local noise induced pattern, except that the latter involves all the antenna's 
oscillation eigenmodes rather than just the monopole/quadrupole harmonics, characteristic 
of generic GWs. This difference proves however insufficient to efficiently filter the stochastic 
GWs from random local disturbances. 

We are thus led to conclude that, in the absence of some a priori knowledge about the 
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GW spectral density Sh(w), the multimode character of a spherical detector does not offer 
a useful alternative to the traditional cross-correlation between two (or more) detectors as 
a recipe to get rid of local noise effects. However, if such a priori knowledge is available 
then the spherical detector naturally yields one order of magnitude better performance (in 
energy) than single readout devices, such as bar and interferometric GW detectors. 
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Appendix A 

A suitable representation for the matrices E^ and E^ m ^ used in equation (|9|) is the 
following: 
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where f2 is a unit vector, and Yi m (Q) are spherical harmonics. 
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